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ABSTRACT 

We caUbrate reverberation-based black hole masses in active galactic nuclei (AGNs) 
by using the correlation between black hole mass, Mbk, and bulge/spheroid stellar 
velocity dispersion, o"*. We use new measurements of cr=K for 6 AGNs and published 
velocity dispersions for 10 others, in conjunction with improved reverberation mapping 
results, to determine the scaling factor required to bring reverberation-based black hole 
masses into agreement with the quiescent galaxy Mbh~o'* relationship. The scatter in 
the AGN BH masses is found to be less than a factor of 3. The current observational 
uncertainties preclude use of the scaling factor to discriminate between broad-line region 
models. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 

The advent of techniques for measuring masses of supermassive black holes (BHs) has led 
to the identification of correlations between the BH mass (Mbh) and various properties of the 
host galaxies. One of the tightest of these relationships is with the velocity dispersion of the 
bulge or spheroid (a*; Ferrarese &: Merritt 2000; Gebhardt et al. 2000a). The objects defining the 
initial MBH~f* relationship were primarily quiescent galaxies, with Mbh determined from stellar 
kinematics or gas dynamics. However, galaxies hosting an active galactic nucleus (AGN), in which 
the BH mass was measured via reverberation mapping (Blandford & McKee 1982; Peterson 1993), 
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have been found to be consistent with foUowing the same correlation (Gebhardt et al. 2000b; 
Ferrarese et al. 2001 [hereafter Paper I]; Onken et al. 2003). 

The values of Mbh derived from reverberation mapping are subject to certain systematic un- 
certainties. In tracing the response of gas in the broad-linc region (BLR) to the variable ionizing 
continuum of the AGN, the time delay between the fluctuations in the continuum and the "rever- 
beration" of the BLR emission lines gives a characteristic radius of the BLR gas, r. The orbital 
velocity at that radius is estimated by the width, AV, of the emission line (specifically, the variable 
part of the line). But the kinematics and geometry of the BLR introduces a scaling factor, /, into 
the reverberation mass equation: 

r 

Mbh = / (1) 

Simple models of BLR morphologies yield / parameters on the order of unity, and the application 
of such estimates, as in the above references, places the previous AGN Mbh-c* data among the 
locus of quiescent galaxies. 

A uniform relationship between the BH mass and properties of the host galaxy on size scales 
beyond the strong gravitational influence of the BH implies a causal connection between the forma- 
tion of the galaxy and the central black hole. Many investigators have explored possible mechanisms 
for how the evolution of the BH and galaxy could be linked (e.g., Silk &: Rees 1998; Hachncit & 
Kauffmann 2000; Adams, Graff, & Richstone 2001; Umemura 2001; Miralda-Escude & Kollmeier 
2003; Merritt Sz Poon 2004), while others have looked for outliers from these relationships as probes 
of the physical drivers of the correlations (e.g. Mathur, Kuraszkiewicz, & Czerny 2001; Wandel 
2002; Bian & Zhao 2004; Grupe & Mathur 2004). 

We present measurements of stellar velocity dispersions for 6 reverberation-mapped AGNs, 
significantly enlarging the sample of objects which can be used to observationally investigate the 
Mbh~<7* relationship for active galaxies. In §2, we describe our observations and analysis method. 
We present our results and discuss the implications for the ensemble average value of the reverber- 
ation mapping scaling factor, (/), in §3. Our conclusions are summarized in §4. 



2. OBSERVATIONS AND DATA REDUCTION 

The velocity dispersions for our sample of AGNs were measured using the near-infrared Ca II 
triplet (CaT). These stellar absorption lines, at rest wavelengths of A8498, 8542, and 8662 A, occur 
in a region of relatively low AGN contribution (Nelson &; Whittle 1995), and are accessible to 
ground-based observations for sources with redshift z < 0.068 (water vapor bands begin to reduce 
atmospheric transparency at longer wavelengths). 

Observations were conducted at Kitt Peak National Observatory (KPNO), Cerro Tololo Inter- 
American Observatory (CTIO), and MDM Observatory. The general observing strategy at each 
of the telescopes was similar. Long-slit spectra of each target were bracketed by quartz-lamp 
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fiat fields and wavelength calibration exposures. The total exposure time for each AGN typically 

exceeded 10,000 seconds. In addition to our targets, we obtained spectra of latc-typc giant stars 
(G8 III-K6 III) to use as spectral templates. Details of the observing runs are given below and in 
Table 1. 

KPNO Observations. We observed at the Mayall 4 m telescope with the Ritchey-Chretien (R- 
C) Spectrograph. The BL380 grating (1200 lines mm~^, blazed at 9000 A) was used with an RG695 
blocking filter and the LBIA CCD. The slit width was 2", with a dispersion of 0.45 A pixel"^ and 
a wavelength range of 8250-9130 A. 

CTIO Observations. We used the V.M. Blanco 4 m telescope with the R-C Spectrograph, 

equipped with the KPGLD-1 grating (790 lines mm" -1, blazed at 8500 A) and an RG665 filter. We 
used the lower right amplifier for the Loral 3K 1 CCD detector within a window of 3072x585 pixels. 
The reduced spectra have a dispersion of 0.83 A pixel"^ over a wavelength range of 7500-10060 A. 
The slit dimensions were 2x344". 

MDM Observations. Our first observations at the MDM 2.4 m telescope were made with 
the MDM Observatory Modular Spectrograph (ModSpec), using a grating of 830.8 lines mm~^, 
blazed at 8465 A. The OG515 order-separating filter was used. The detector for these observations 
was "Charlotte", a thinned, backside- illuminated, SITe 1024x1024 CCD. The spectra covered a 
range of ~1400 A, with a dispersion of about 1.43 A pixel"^, and a slit width of approximately 
2". The second MDM observing run used the ModSpcc setup described above, but utilized a 
different detector. "Echelle", a thinned, backside-illuminated, SITe 2048x2048 CCD was used with 
a windowing of 330x2048 pixels. 

2.1. Data Reduction and Analysis 

We used IRAF^ and XVista^ for the reduction of different subsets of the data. However, the 
basic strategy was the same: the spectra were fiat-fielded, sky-subtracted, and placed on a linear 
wavelength scale. 

The AGN spectra in our sample were often contaminated by the broad O I A8446 emission 
line that appears blueward of the CaT, and these lines were removed with a high order polynomial 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation 

^See http:/ /ganymede. nmsu.edu/holtz / xvista 



-4- 



in order to isolate the CaT absorption lines. None of our objects show evidence of emission from 
high-order Paschen emission lines, which can also contaminate the CaT lines. 

3. RESULTS AND DISCUSSION 

3.1. Measuring Velocity Dispersions 

We measured the velocity dispersions from the reduced spectra with the Fourier correlation 
quotient method (FCQ; Bender 1990). Cross-correlating the galaxy spectra with those of the 
template stars produces values for the central velocity dispersions. 

Our galaxy spectra and the best template-star fits from the FCQ routine are shown in Figure 1. 
The results of our velocity dispersion analysis are listed in Table 2. As noted in Paper I, the formal 
errors reported by the FCQ fitting routine tend to underestimate the total uncertainties for the cr* 
measurements. Thus, in cases where the FCQ output errors were smaller than 15%, the error bars 

have been brought up to this threshold. 

Based on visual inspection, the least satisfactory fit from Figure 1 is that of Akn 120. However, 
consistent (within 1 a) results are found from three other methods: (1) a simple Gaussian fit to 
the strongest absorption line; (2) the second moment of the profile for that line; and (3) an IDL 
routine for penalized pixel fitting^. 

We examined whether the AGN continuum could dilute the stellar absorption lines in a manner 
that would affect the a* measurements. This was tested by artificially including an additional 
continuum contribution and running the FCQ routine again. Over a large range of continuum 
levels, no significant change was seen in the resulting cr* values. 

Additional stellar velocity dispersion data for galaxies hosting AGNs have been taken from the 
literature. 

3.2. Measuring Virial Products 

New analysis of reverberation mapping data has produced updated measurements of the virial 
products (r AV'^/G) for a sample of 35 AGNs, including the 16 which now have velocity dispersions 
(Peterson et al. 2004). These virial products, which, from equation (1), can also be written as 
Mbh//, differ from the majority of previous reverberation mapping analyses in the calculation of 
AV. Whereas earlier work typically characterized the velocity width via the full-width at half- 
maximum (FWHM) of the emission line, Peterson et al. (2004) have found more consistent virial 
products when using the line dispersion (i.e., the second moment of the line profile; crime) instead. 



^See Cappellari & Emsellem (2004), and http://www.strw.leidenuniv.nl/~mcappell/idl/ for details. 
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This will have important implications for comparing our results with previous publications. 



We assume that the Mbh-ct* correlation found in quiescent galaxies also holds for AGNs and 
their host galaxies. We can then determine the factor required to scale the AGN virial products 
(Mbh//) to the quiescent galaxy relation (holding cj* fixed at the measured value); that is, we 
calculate the average scaling factor, (/), that, when multiplied by the virial products, brings the 
AGN Mbh^o"* relation into agreement with the quiescent galaxy relationship. 

Because published determinations of the slope of the quiescent galaxy relation have yet to 
converge, we selected the most prominent slope values from either end of the quoted range. Thus, 
we consider both the Tremaine et al. (2002; hereafter T02) slope of 4.02, and the Ferrarese (2002; 
F02, henceforth) slope of 4.58, and determine separate scaling factors, (/t) and (/p), respectively. 

While questions remain as to the reliability of the current fits to the quiescent galaxy Mbh-c* 
relation, especially with respect to the accuracy of the black hole masses from stellar kinematics 
(see Valluri, Merritt, & Emsellem 2004; Cretton k. EmscUcm 2004; Richstone et al. 2004), we hope 
our use of the two slope values will give a more reliable estimate of the uncertainties involved while 
also avoiding entanglement in the continuing controversy. 

The best fit to the AGN Mbh^^* values was accomplished with the orthogonal regression 
program, GaussFit^ (version 3.53; Jcffcrys, Fitzpatrick, &: McArthur 1988), which accounts for 
errors in both (Mbh//) and a^,. For the purposes of our GaussFit analysis, the asymmetric variances 
in the virial products were symmetrized as the mean of the upper and lower variances. A fit was 
made to the equation 



where a is the normalization of the relationship, (3 is the (fixed) value of the slope, and where 
(Mbh//) is in solar masses and a* is in km s~^. For each slope, the best-fit value of a was 
calculated. The which was minimized to determine the best fit is given by 



where ayii and a^i are the uncertainties in the virial product and the velocity dispersion, respec- 
tively, for the ith data pair (Press et al. 1992). The virial products and velocity dispersions used 
in our fits are given in Table 3. 



3.3. The Mbh-o-* Relation 




(2) 




(3) 



''GaussFit is available at ftp://clyde.as.utexas.edu/pub/gaussfit/ 
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The published values of the quiescent galaxy intercept correspond to ap = 8.22±0.08 and 
aT = 8.13ib0.06. Our fits with fixed slopes are shown in Figure 2, and yield ct^*^^ = 7.48ib0.13 
and a^^^ = 7.39=b0.12. Remarkably, both slopes give the same best fit value, log(/) = 0.74. 
Hence, (/f) = 5.5±1.9 and (/t) = 5.5±1.7. The psr degree of freedom, xt^ for these two fits 
are 2.90 and 2.87, respectively. 

Table 4 shows the resulting a and /? values for a variety of fitting constraints. When we allow 
P to vary in the fit to the AGNs, the result is consistent with both the T02 and F02 slopes, albeit 
with a large uncertainty. The distribution of points with previous cr* measurements in Figure 2 
differs from that of Paper I because of the improved virial product data from Peterson et al. (2004) . 

Because the reverberation mapping results for IC 4329A and NGC 4593 provide only upper 
limits on the black hole masses, these data were omitted from the above analysis. The formal error 
bars on Mbh can extend below zero because the uncertainty of the time delay is not restricted to 
positive values. If we take the uncertainty in log Mbh to be aM/{MlnlO), we can avoid taking 
the logarithms of negative numbers, and these two data points can be included in the fit. The 
uncertainties are large enough, however, that the fit is not changed in a statistically significant way. 
Additionally, the fit is not strongly affected by the exclusion of NGC 4051 (see Table 4), which 
seems to be an outlier in the AGN radius-luminosity relationship (see the Appendix of Vestergaard 
2002). When we do not include NGC 4051, we find (/f)=5.1 and (/t)=5.2, differences of less than 
8% from our best fit value of 5.5. 

The most common assumption for converting a virial product to Mbh is that of Netzer (1990) 

Mbh = (4) 

which implicitly assumes an isotropic velocity distribution and AV = Vfwhm / 2, where Vfwhm is 
the line full width at half maximum. Peterson et al. (2004) measure the second moment of the line 
profile, (Tiine, rather than VpwHMi so an isotropic velocity field simply gives / = 3 with no other 
assumptions. Denoting the scale factor as e when using VpwHM and / when using cxiinc, we find 
that whereas Netzer's assumptions give e = 0.75 (as in eq. [4]), the value of (/) we derive above 
implies (e) = 1.4, i.e., black hole masses ~1.8 times higher^. Our masses are only 62% as large as 
those expected from the criine-VpwHM relation found in a model-dependent fit to the distribution 
of Vfwhm in a sample of quasars and Seyfert galaxies (McLure &; Dunlop 2001). 

Applying our derived value of (/) = 5.5 (or (e) = 1.4) will remove the systematic bias between 
the virial product and mass of the black hole (modulo the accuracy of the particular Mbh-ct* fit 

chosen). Yet, estimates for individual AGN black hole masses may fall substantially off the i\fBH~c* 
relation, as there are still rms deviations from the fits in the mass direction of factors of 2.9 (F02 

slope) and 2.6 (T02 slope). 



®Our value of e is calculated assuming crune = Vfwhm /2; for typical line profiles, changing the shape of the line 
alters the crune-VFWHM conversion by an amount smaller than the observed scatter in the Mbh-ct, data. 
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Whether the scatter in Mbh is due to differences in the BLR properties in each AGN (cor- 
responding to slightly different / values for each object) or a problem in our assumption that all 
AGNs fall on the quiescent galaxy relationship remains unclear. However, models of the BLR may 
shed light on this question. 



3.4. BLR Models 

Given our determination of (/), we can try to learn something about the structure and kine- 
matics of the BLR. Modeling of the BLR, taking into account different radial profiles and discrete 
time sampling of example light curves, will be undertaken in a later paper. Here we describe only a 
thin, rotating ring. In this case, / can be thought of as relating the Keplerian velocity (K-ot) at the 
radius of the ring (which is traced by the appropriate time delay) to the observed line dispersion: 

f=(^)\ (5) 

For a ring of inclination i (where i=0° for face-on), this model gives / = 21n2 / sin^z. 

Working within a thin-ring model context, Wu & Han (2001) calculated the inclinations neces- 
sary to scale published reverberation mapping results for 11 AGNs to the quiescent galaxy Mbh-ct* 
relation. With our revised virial product results, we find that this very simple model fails for two 
of our 16 AGNs. NGC 3783 and Mrk 110 would require sinz values greater than unity to scale 
them to the quiescent galaxy Mbh^c* relation. For our other sources, we do not find any evidence 
to support the inclination trends with line width or radio loudness that were claimed by Wu & 
Han. Further constraints on the BLR geometry are not feasible given the sizes of the current 
uncertainties. 



3.5. Gravitational Redshifts 

In principle, an independent measure of the AGN black hole mass can be obtained by de- 
tection of gravitational redshifting of the emission lines (e.g., Netzer 1977; Peterson et al. 1984). 

Kollatschny (2003b) finds that the variable parts of the strong broad lines in the optical spectrum 
of Mrk 110 are redshifted relative to the systemic velocity. The reverberation time lag and redshift 
of each line relative to systemic are anticorrelated and appear to be consistent with a gravitational 
redshift caused by a central mass Mgrav ~ 1-4 x 10^ Mq, rather higher than our reverberation-based 
mass^ (Table 3) of Mbh = 2.5 x 10^ Mq. 



^Kollatschny (2003a) obtains a reverberation-based mass estimate of 1.8 x 10^ Mq, but using FWHM as the line- 
width measure and e = 1.5 (see §3.1). The FWHM data of Peterson et al. (2004) gives a result consistent with that 
of Kollatschny, however our preferred method of measuring crune and using the empirical calibration of {/) yields a 
black hole mass approximately 40% larger. 
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Regardless of the particular value of / appropriate to Mrk 110, the offset between the reverberation- 
based virial product and the gravitational redshift mass is in the direction expected. However, the 
value of Mgrav places Mrk 110 even further above the Mbh-c* relation, which predicts Mbh ~ 4 x 
IO^Mq. 



If the BLR is modeled as a simple disk, then 

(Mbh//) sin^i 



Mgrav 2 In 2' 



(6) 



and the difference between these two measurements allows us to infer that i ~ 30°, slightly larger 
than the 21 ± 5 degrees found by Kollatschny (2003b). 

Measuring gravitational redshifts appears to be a promising technique for independently es- 
timating central masses. However, applying the method is non-trivial because the gravitational 
redshifts are small, typically only several tens to a few hundred kilometers per second. Moreover, it 
is not a shortcut around reverberation mapping because (a) the redshift must be measured in the 
variable part of the emission line (i.e., the root-mean-square spectrum formed from the monitoring 
data; see Peterson et al. 2004) and (b) the time lag, which can vary with time, needs to be measured 
simultaneously. The spectral resolution required to measure such small redshifts reliably is higher 
than usually employed in reverberation mapping campaigns, with the exception of Kollatschny's 
(2003a) program on Mrk 110 and the Hubble Space Telescope monitoring program on NGC 5548 
(Korista et al. 1995). A re-examination of the NGC 5548 spectra from the latter campaign shows 
some evidence for redshifts at the expected levels, although the errors are quite large. Additional 
data will thus be required to examine this method of black hole mass measurement more thoroughly 
and to determine whether Mrk 110 is truly a significant outlier from the Mbh~<7* relation. 



3.6. 0-* Versus FWHM([0 III]A5007 A) 

Several recent studies have suggested the use of the FWHM of the [O III] A5007 A emission line 
as a proxy for a*. Nelson &; Whittle (1996) examined a sample of 66 Seyfert galaxies with both a* 
and FWHM([0 III]) measurements and found scatter of 0.20 dex around a 1:1 correspondence (see 
their Figures 5 and 7a). Based on these trends. Nelson (2000) and Shields et al. (2003) used [O III] 
data as a substitute for a^, in AGNs out to z~3. Boroson (2003) applied the same arguments to data 
from the Sloan Digital Sky Survey Early Data Release and found that using FWHM([0 III]) could 
reproduce an MBH^O'*-like relationship, although with a larger scatter than has been found for 
quiescent galaxies (i.e., a factor of ~5). Statistically identical distributions in FWHM([0 III]) were 
found for the broad-line and narrow-line Seyfert Is in an X-ray selected sample of AGNs, which 
placed the narrow-line Seyfert Is preferentially below the Mbh-c* relation (Grupe &; Mathur 2004). 

With the FWHM([0 III]) data tabulated by Nelson (2000) for our 16 AGNs, we looked at 
the relationship with cr*. Figure 3 shows that 25% of the sources have FWHM([0 III]) data that 
deviate by >0.2 dex from the values expected based on their velocity dispersions. There is no 
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evidence of the correlation between the discrepant objects and their radio power that was found 
for a larger sample of Seyferts (Nelson & Whittle 1996), indicating that these differences may not 
be due to a systematic acceleration of the [O IIl]-emitting gas by the radio source. Overall, the 
weighted mean difference from equality is 0.03 dex (in the sense of larger FWHM relative to a^) 
with an rms scatter of 0.15 dex. In the way of commentary, we simply echo sentiments expressed 
by others that while use of such a proxy may be valid for a large sample of AGNs, it can also fail 
dramatically for individual objects. 

4. CONCLUSION 

With the addition of 6 new velocity dispersion measurements for rcvcrbcration-mapped AGNs 
and making use of improved reverberation mapping results, we tie together the Mbh^o"^ relation- 
ships for quiescent galaxies and AGNs. This allows us to calculate the average scaling factor, (/), 
which removes the statistical bias between the virial product generated by reverberation mapping 
and the black hole mass. For the F02 and T02 fits to the quiescent galaxy Mbh-c* relationship, 
we find (/) = 5.5ibl.9 and 5.5ibl.7, respectively. These values of (/) apply specifically to virial 
products using the dispersion of the emission lines, rather than measurements of the FWHM. While 
modeling of the BLR and studies of emission line gravitational redshifts may eventually lead to a 
better understanding of the structure and kinematics of the BLR, further work is needed. 
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Fig. 1. — Normalized spectra of the CaT region for the 6 AGNs in this study. Each spectrum has 
been offset in flux for clarity. The dashed line indicates the best fit obtained with the FCQ method. 
The spectrum of Akn 120 is truncated where the FCQ fit becomes unstable. The spectra shown 
for Mrk 590 and NGC 7469 are from KPNO; IC 4329A and NGC 3783 are from MDM. 
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Fig. 2. — Virial product, Msu/f, versus host galaxy velocity dispersion, o"*. Solid points indicate 
objects with measurements presented here. Open points represent AGNs with previously pub- 
lished C7* data. The solid line indicates the F02 slope of 4.58, with the y-intercept shifted downward 
by (/f)- The dotted line denotes the T02 slope of 4.02, with the y-intercept scaled down by (/t)- 
The vertical scale on the right uses our derived offset: (/f)=(/t)=5.5. 




Fig. 3. — cr* versus FWHM([0 III]) for our 16 AGNs. The solid line shows a 1:1 correspondence 
between cr* and the equivalent for a Gaussian line profile, FWHM([0 III]) / 2.35. 
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Table 1. Observing Log 











Resolution 


Targets 


Run 


Telescope 


Instrument 


UT Dates 


(km s"^) 


Observed 


MDM-a 


MDM 2.4m 


ModSpec 


2001 Oct 20-24 


95 


Mrk 590, NGC 3227, NGC 7469 


KPNO 


KPNO 4m 


R-C Spec 


2001 Oct 29-Nov 1 


60 


Akn 120, Mrk 590, NGC 7469 


MDM-b 


MDM 2.4m 


ModSpcc 


2003 Mar 12-13 


95 


IC 4329A, NGC 3783 


CTIO 


CTIO 4m 


R-C Spec 


2003 Apr 18 


90 


IC 4329A, NGC 3783 



Table 2. Velocity Dispersion Measurements 







Redshift 






Galaxy 


R(-<lsliift 


Rcfort'iiccs 


Run 


(km s ^) 


Akn 120 


0.032296±0.000143 


1 


KPNO 


239±36 


IC 4329A 


0.016054±0.000050 


2 


CTIO 
MDM-b 


131+2° 
121±18 


Mrk 590 


0.026385±0.000040 


3 


KPNO 

MDM-a 


201±30 

188±28 


NGC 3227 


0.003859±0.000010 


3 


MDM-a 


139±21 


NGC 3783 


0.009730±0.000007 


4 


CTIO 

MDM-b 


87±13 

108±16 


NGC 7469 


0.016317±0.000007 


5 


KPNO 
MDM-a 


149±22 
157±24 



References. — (1) Falco et al. 1999; (2) Willmer et al. 1991; (3) de 
Vaucouleurs et al. 1991; (4) Theureau et al. 1998; (5) Keel 1996. 
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Table 3. Mbh-ct* Data 





Virial Product'' 


Black Hole Mass'' 


cr. (avg) 




GalEixy 


(10« Mq) 


(10« Mq) 


(km s"^) 


References 


3C 120 


lO.lti;! 




162±24 


2 


3C 390.3 


52.2±11.7 


289±64 


240±36 


3 


Akn 120 


27.2±3.5 


150±19 


239±36 


1 




1 sn+^-^^ 

i.»U„2.16 




1 994-1 Q 


1 

L 


Mrk 79 


9.52±2.61 


52.4±14.4 


130±20 


4 


Mrk 110 


4.57±1.10 


25.1±6.1 


86±13 


4 


Mrk 590 


8.64±1.34 


47.5±7.4 


194±20 


1 


Mrk 817 


8.98±1.40 


49.4±7.7 


142±21 


4 


NGC 3227 


7.67±3.90 


42.2±21.5 


131±11 


1,2 


NGC 3516 


7.76±2.65 


42.7±14.6 


164±35 


5 


NGC 3783 


5.42±0.99 


29.8±5.4 


95±10 


1 


NGC 4051 


0.348±0.142 


1.91±0.78 


84±9 


2,4 


NGC 4151 


2.42±0.83 


13.3±4.6 


93±14 


4 


NGC 4593^^ 


0.9811;™ 


c 0(7+9. 37 


124±29 


2 


NGC 5548 


12.20±0.47 


67.1±2.6 


183±27 


4 


NGC 7469 


2.21±0.25 


12.2±1.4 


152±16 


1 



''Prom Peterson et al. (2004). 
''Scaled using / = 5.5. 
''Excluded from fits. 

References. — (1) This work; (2) Nelson & Whittle 1995; (3) Green et al. 
2003; (4) Paper I; (5) Arribas et al. 1997. 



Table 4. Mbh-o"* Fitting Results 



Constraint 


Slope (/3) 


Intercept (a) 


xl 


(/> 


F02 slope 


4.58 


7.48±0.13 


2.90 


5.5±1.9 


T02 slope 


4.02 


7.39±0.12 


2.87 


5.5±1.7 


none 


4.11±1.07 


7.40±0.21 


3.11 


N/A 


F02 slope, no NGC 4051 


4.58 


7.50±0.14 


3.09 


5.2±1.9 


T02 slope, no NGC 4051 


4.02 


7.42±0.12 


2.96 


5.1±1.6 



